Simplifying Kaufman's Solution of the Two-Dimensional Ising Model 
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We considerably simplify Kaufman's solution of the two-dimensional Ising model by introducing 
two commuting representations of the complex rotation group SO(2n,C). All eigenvalues of the 
transfer matrix and therefore the partition function are found in a straightforward way. 
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Since Onsager's solution Q in the transfer matrix ap- 
proach g] of the two-dimensional Ising model jjj with 
vanishing magnetic field and its subsequent simplifica- 
tion by Kaufman ||, there have been a number of re- 
lated as well as alternative solutions, see e.g. Baxter's 
book on exactly solved models in statistical mechanics 
|^| and references therein. Among the transfer matrix 
solutions are the ones by Schultz, Mattis and Lieb, by 
Thompson, by Baxter and by Stephen and Mittag ||. 
Nevertheless, the author of this work feels that there is 
still room for a nice and straightforward solution. A com- 
pletely self-contained and detailed account of this work 
may be found in 0. 

We study the two-dimensional Ising model with zero 
magnetic field on a square lattice with m rows and n 
columns subject to toroidal boundary conditions. The 
transfer matrix is expressed in terms of the generators of 
two commuting representations of the complex rotation 
group SO(2n, C). These representations naturally arise 
from projected bilinears of 2™ x 2™ spin matrices. Con- 
servatively speaking, we reduce Kaufman's approach to 
its essential steps, avoiding in particular the doubling of 
the number of eigenvalues of the transfer matrix and sub- 
sequent rather involved arguments for the choice of the 
correct ones. Additionally, there is no need to investigate 
the transformation properties of the spin matrices. 

Our notation is in the spirit of || , with sans serif cap- 
itals reserved for 2" x 2" matrices. The structure of this 
work is as follows: After defining the model and its trans- 
fer matrix T, we express T in terms of 2™ x 2™ spin ma- 
trices X„, Y u , Z„. A rescaled transfer matrix V is de- 
fined whose eigenvalues are, up to a trivial factor, the 
eigenvalues of T. We define further spin matrices V v and 
two commuting projection classes Jjtg and of their 
bilinears. After investigating the relevant properties of 
the J^o, we express V in terms of them. We introduce 
2n x 2n matrices J a p whose algebra is identical to that of 
Jj^g and define matrices in terms of the J a( g such that 
the relation between and J a p is closely related to that 



between V and J~j. The result of the well-known diago- 
nalization procedure for the is given and the analogy 
between and V exploited for the diagonalization of 
V. The eigenvalues of V and the partition function are 
found explicitly. 

The energy is given by E = E a + Eb with 

m n m n 

E a = J a J]] J]] S/ii/S^+l,!/, Eb = Jb J]] J]] SfiuS^u+l, 
/i— 1 u— 1 M— 1 f=l 

with j3 1 = kT . J a and Jb are temperature-independent 
interaction energy parameters that we assume to be neg- 
ative. We identify rows 1 and m+1 and columns 1 and 
n+1, i.e., the lattice is wrapped on a torus. The s M „ can 
take the values ±1. The partition function is then given 
by 



Z(a,6) = ^--.^exp(-/3£), 



(2) 



with the definitions a — —[3J a and b — —flJb- Z can be 
expressed with the help of a 2™ x 2™ transfer matrix T, 
Z(a,b) = TrT m with T defined by its elements (s n+ i = 
si), 



(7r|T|7r') = exp(as v s' v + bs v s v+1 ), 



(3) 



u=l 



where 7r M = . . . , s M „} for /i = 1, . . . , to. We can split 
T into a product of two matrices T = V^V^, defining 
and Vf, by their elements 



n 

<7r|V b |7r'> = Y[ exp(&s„s v+ i). (5) 



u=l 



With the help of the Pauli matrices a x , a y , a z , and the 
2x2 unit matrix 1, we define Hermitian 2™ x 2™ spin 
matrices by the direct products 
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and analogously for Y v and in terms of a v and a z , 
respectively. With a > denned by sinh(2a) sinh(2a) = 
1, we can write = [2 sinh(2a)]™ /2 V a with 



and 



V ( , = Y[ CX P (aX„) 



Vb = II cxp (bZvZv+i) . 



(7) 



(8) 



where we have identified 7- n +\ = Z\. The transfer matrix 
may then be expressed as T = [2 sinh(2a)]"/ 2 VftV a . Due 
to the cyclic property of the trace, we may rewrite the 
partition function (|^) as 

Z(a,b) = [2sinh(2a)] m "/ 2 TrV m , (9) 
where V is defined by the Hermitian matrix 



with 



V = V a/2 V b V a/2 



V a/2 = n exp(aX,/2) 



(10) 



(11) 



so that V 2 ^ 2 = V a . If Afc are the 2 n eigenvalues of V, we 
have 



Z( a ,6) = [2sinh(2 a )]" m / 2 ^A1». 



(12) 



k=l 



Our task is therefore to find the eigenvalues of V. 



which anticommutes with every T^, {T^, Ux} = 0. Now 
we can write for the matrices appearing in the exponents 
of Eqs. (!) and (111]) 



X,, = --[^zu, r 2w -i], v = l,...,n, (17) 
Z V Z V+1 = —-[V 2v +i,r 2u ], v=l,...,n-l, (18) 



z„z x = -Ux[ri,r an ]. 



(19) 



So far our treatment has been rather similar to Huang's 
write up Q of Kaufman's approach . Our subsequent 



treatment rests on the observation that the formulation 
@ of V with Eqs. (|) and (|TTJ) involves only the product 
Ux of all r„ and bilinears r a r<3, see Eqs. This 
will allow us to express V in terms of the elements of two 
commuting algebras of projected bilinears of the I" „. 
With the help of the projectors 



p ± = i(l±u x ), 



define the matrices 



so that 



Since J 



Jq/3 — J Q/3 + J a/ 3, 



J ± =±J ± 



(20) 



(21) 



(22) 



Define the 2n matrices (y = 


l,...,n) 




v a/2 


^1v-\ — Xi ■ • 


' Xy_lZy, 


(13) 




^2v = Xi • ■ 


• X^-iY^, 


(14) 


with 


which obey 








{r M ,r,} = 




(15) 


and 


Define further the matrix 






U x = X 1 .-.X n = fT 1 IV- 


• t~2n, U x = 1, 


(16) 





j^g = — J^ Q , there are n(2n — 1) such independent 

matrices of each kind J^g and J"^. It is straightforward 
to show that their algebra decomposes into two commut- 
ing parts, [J„0,J~$] = 0, which obey identical algebras 

^ = i( s ^% + ^ J «7 - - S 0J £ S ). (23) 

Next note that with Eqs. @, and @ we 

can write 

X, = 2(J+^_ 1 +J 2 -, 2 „_ 1 ), i/ = l,...,n, (24) 
Z„Z„ +1 = 2(J+ +1 2 ,+J 2 - +1 2 J, i/ = 1, . . . , n-l,(25) 
Z n Zi = -2U X (J^ 2 „+Jr,2n) = -2(JU-Jr,2„)- (26) 

This allows us to express V a / 2 from Eq. ( |lT| ) and V(, from 
Eq. (||) in terms of the J^g, 

n 

V Q/2 = 1] exp[a(Ji l2v -i + Ji,,a„-i)] = V+ /2 V a - /2 (27) 



V ^/2 = J! eX P( aJ 2^,2^_ 1 ) = 



V fc - exp[-26(J+ 2n -Jr 2 J 



(28) 



JJ exp[26(J+ +lj2iy +J 2 



= v+v- 



with 



n- 1 



V± = cxp( T 2feJ± 2 J J] exp(26J± +lj2 J. 



(29) 



(30) 



i/=i 



The rescaled transfer matrix V defined in Eq. (|10| ) reads 
then V = V+V" with 



V* — V+ V+V+ 

V a/2 V 6 V a/2' 



[v+,v-] = o. 



(31) 
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Define N x N matrices J a p by their elements 



(32) 



where Greek and Latin indices run from 1 to N. Since 
J a p = —Jg a , there are N(N — l)/2 such independent 
matrices. As can be easily checked, they obey the algebra 
(p3|), if we set N = 2n. Now consider the matrices S = 
exp(icap Jap), where c a g are arbitrary complex numbers. 
The matrices S form the group SO (TV, C) of complex N x 
N matrices with = det S = 1. 
Define the SO(2n, C) matrices 



v ± = v a/2 v±v a/2 



with 



V, 



a/2 



'. [ exp(aJ 2l ,.2i^-i) 



and 



n-l 



= exp(=F26Ji :2 «) Y\_ cxp(2bJ 2l s+i,2v), 



(33) 



(34) 



(35) 



u=l 



in analogy with V ± , V± /2 , and in Eqs. (|l|), (||), and 

(|30|). Since a and 6 are real, the matrices V a /2, V^, and 
are not only orthogonal, but also Hermitian, so the 
V r± have only real eigenvalues and in each case a complete 
set of orthonormal eigenvectors. 

It is well known ||, |j how to diagonalize matrices 
of the types V ± . Applying similarity transformations 
V^ 1 = S±V :iz S± 1 with certain explicitly known matrices 
S±, one obtains 



V s =exp 



with 7^ defined by 




(36) 



cosh jk — cosh 2a cosh 2b — cos — sinh 2a sinh 2fe. (37) 

We fix the sign of 7& by defining 7^ = 2a for 6 = and 
then analytically continuing to other values of b. For 
k = 1, . . . , 2n — 1, this means 7^ > 0. On the other hand, 
for 70 this means 



7o = 2(5 - b). 



(38) 



Our sign convention for the jk and in particular for 70 
allows us to treat all jk on an equal footing for both a > b 
and a < b, i.e., irrespective of the temperature. Note 
that the are only block diagonal with 2x2 blocks. 
It would be trivial to diagonalize Vg , but the form given 
by Eq. (^6|) is most convenient for our purposes. 

It is straightforward to show that the matrices S± are 
elements of SO (AT, C) and may therefore be written as 



S± = exp{ic a pJ a p). Now use the same parameters c a/3 
to define the 2™ x 2 ra -dimensional transformation matrix 

S = S+S_, S± = exp^c^J^,), (39) 

and write 

V s = SVS" 1 = S+V+S^S-V-SI 1 = V+V s . (40) 

The factors defining Vg in terms of a, b, c^, and the 

have the same structure as the in terms a, b, c^g, 
and the J a p- Now imagine using the Baker-Campbell- 
Hausdorff formula || 

exp(A) exp(i?) = 

exp(A+B-±[B, A]+±{[A, [A, B]] + [B, [B, A]]} + • • •) 

(41) 

to work out all products of exponentials in Vg . Since the 
JLg, Jig, and J a g obey identical algebras, the result is 



so that 
V s 



= exp 



ex P Eto J a^-i ' 



(42) 



exp 



v=l 

1 - 

t^>-i(1 + Ux)X„ 



4 

v=l 



n 

Z E^-2(1-U X )X, 



(43) 



with the same 7^ as defined in Eq. (|37|) and the subse- 
quent sign convention. 

To diagonalize Vs, define another similarity transfor- 
mation Vy = RyVsRy 1 with Ry and its inverse given by 



R± i =2 -"/ 2 n(i±*Y,). 



(44) 



Since RyX^R,, 



Z„, this transformation takes Vs into 



Vy = RyVsRy 
1 , 

= exp 



1 " 1 " 

j E 72,-i(l+Uz)Z„+ z £ 72,- 2 (l-Uz)Z, 

(45) 



with Uz = Zi • • • Z„. 

The matrix Vy is diagonal, but we still have to deter- 
mine its elements. Uz is a diagonal matrix with elements 
+1 and —1 occurring in equal numbers. For each element 
holds, if an even (odd) number of Z„ provides a factor 
— 1, the matrix element of Uz is +1 (— 1). This means (i) 
a matrix element of (1 + Uz)/2 is 1 (0) if an even (odd) 
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number of 7- v provides a factor — 1, (ii) a matrix element 
of (1 — Uz)/2 is 1 (0) if an odd (even) number of 7- v pro- 
vides a factor —1. It follows that the 2™ eigenvalues of V 
split into 2" -1 eigenvalues of the form 

exp^^(±)72„-ij ' (46) 
and 2 n_1 eigenvalues of the form 




where in the first (second) case all sign combinations with 
an even (odd) number of minus signs occur. This is re- 
flected by the indices "e" and "o" in our result for the 
partition function, 



Z(a,b) = [2sinh(2a)] m ™/ 2 
1 



5Z ex P T ^2(±)l2v-i +^exp y ^(±) 72l ,_ 2 



[2 sinh(2a)] 



rnn/2 



n n n n 

n [2cosh(^ 72fe -i)] +n [2sinh(| 72fe -i)] +n [2cosh(| 72fe - 2 )] - n [ 



fe=l 



fc=l 



fe=l 



2 sinh ( — 72fc- 2 



(48) 



The last term within the braces has a sign differing from 
that in This is due to our different sign conven- 
tion for 7o. The eigenvalues of T are of course obtained 
by multiplying Eqs. ( [46| ) and (]47]) with the trivial factor 
[2sinh(2a)]™/ 2 . 

The results for the eigenvalues of T and the partition 



function are the starting point for the analysis of the 
thermodynamic properties of the two-dimensional Ising 
model, the most interesting case being the thermody- 
namic limit m, n — > oo. Such analyses can now proceed 
as usual (see, e.g., H, Q, ||) and will not be repeated here. 
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